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■ Abstract. The launching of up to four astrometric missions in the next 

fNJ , decade will enable us to make measurements of stars in tidal streamers 

from Galactic satellites of sufficient accuracy to place strong constraints 
on the mass distribution in the Milky Way. In this paper we simulate 
observations of debris populations in order to assess the required prop- 
erties of any data set chosen to implement this experiment. We apply 
our results to find the desired target properties of stars (e.g., accuracy 
of velocity and proper motion measurements) associated with the dwarf 
spheroidal satellites of the Milky Way. 
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O 1. Introduction 

o ' 

Stars lost from a satellite galaxy in orbit around the Milky Way will gradually 
drift ahead or behind their parent system to form streams of debris along its 
Q-f orbit (Tremaine 1993; Johnston 1998). Assuming we know the phase-space co- 

ordinates of the parent satellite, the unique properties of tidal debris streams 
can be exploited to measure the potential of the Milky Way in two ways: 
MD Algorithm: Integrate the satellite's orbit backwards and forwards and 
try to fit any available observations to the orbit, with the "best" potential cor- 
responding to the closest fit (e.g., Murali & Dubinski 1999). 
JZSH Algorithm: Integrate the stars and satellite backwards and choose the 
"best" potential to be the one in which most stars recombine with the satellite 
(e.g., Johnston et al. 1999, hereafter JZSH99). 

The MD algorithm can be applied with limited observations (e.g., just angular 
positions and line-of-sight velocities) of any accuracy. However, in reality debris 
in trailing/leading streamers are not on exactly the same orbit as their parent 
satellite, but rather in orbits offset in energy from the satellite's orbit by ~ ±e 
where 

e = r tide-^: « sv civc . (1) 

Here, rtide = rs is the tidal radius of the satellite, ^(r) is the Galactic potential, 
■Ucirc is the circular velocity of the Milky Way and s is the tidal scale 
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where m sa t is the satellite's mass and Mo a i is the mass of the Galaxy enclosed 
within the satellite's orbit. Hence, once the observational errors are small enough 
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to detect the offset of the particles from the satellite's orbit the MD algorithm 
will be biased by this systematic offset, rather than limited by the errors them- 
selves. In contrast, the JZSH algorithm exploits the fact that debris stars are 
on orbits with different time periods than the satellite's own and hence will only 
work with observations of five of the phase-space coordinates of sufficient ac- 
curacy that the offset from the orbit can be detected (the sixth coordinate, if 
poorly known, can always be found by applying the principle of energy conser- 
vation). The dividing line between useful application of these two approaches 
will therefore occur where the velocity measurement errors Av satisfy 

v c i rc Av ps e => ps s. (3) 

^circ 

Once measurements of proper motions and velocities of debris can be made 
with accuracies of order 3-30 km/s the second approach will be more accurate 
for all Milky Way satellites (depending on the mass and orbit of the satellite 
- see Table 1) Such accuracies will become possible with the next generation 
astrometric missions: NASA's Space Interferometry Mission (SIM) and ESA's 
Global Astrometric Imager for Astrophysics (GAIA). 

JZSH99 have already demonstrated that the JZSH algorithm can in prin- 
ciple recover both the mass and geometry of the Milky Way with few percent 
accuracies using just 100 stars. In this paper, we examine how these uncertain- 
ties grow with the limitations likely to be present in a real data set. 

2. Methods 

We test how well we can recover the Milky Way potential by "observing" debris 
populations at the end of simulations of satellite destruction and then applying 
the JZSH algorithm to our data sets. We illustrate our results in this paper with 
tests on simulations referred to as Models 1/2/3/4, which had satellite masses 
5.1e6/2.9e7/6.9e7/1.0e8 M , on orbits with pericenters at 43/43/26/30 kpc and 
radial time periods of 2.5/2.5/2.0/1.0 Gyears. These parameters corresponded 
to tidal scales s = 0.026/0.046/0.075/0.081 (see equation [3]). The simulation 
technique and Galactic model are described in Johnston, Spergel k, Hernquist 
(1995). 

We restrict ourselves to asking how well we can recover the circular velocity 
of the halo component that was used in the simulations, assuming we know the 
other parameters in the potential. As an an indication of the error, we attempt 
to recover v c i rc from ten data sets, each containing a different group of particles 
or a different error realization. We take the "error" to be the dispersion in the 
recovered i> c i rc in units of the true circular velocity (= o- Vcirc /v c i TC ). 

3. Results 

Measurement Errors of Debris Stars As noted in §1, we expect the JZSH 
algorithm to require knowledge of velocities such that Av/v c i rc < s. To test this, 
we observed 128 stars in Models 1-4. We then added errors chosen from gaussians 
of mean zero and dispersions Av and Av/d to each line of sight velocity and each 
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Figure 1. Fractional error in f c j rc recovered as a function of velocity 
measurement errors (left hand panel) and number of orbits experienced 
by debris star (right hand panels) 



component of proper motion respectively, and ran the JZSH algorithm to recover 
w c i rc . This was repeated 10 times with the same sample of stars (but different 
error realizations) for Av = 1, 5, 10, 15, 20, 25, 30km/s and the dispersion among 
the results from the 10 realizations calculated. The errors found, shown in the 
left hand panel of Figure 1, confirm the intuition outlined in §1. 

Length of Tidal Streamers A typical debris star will drift by 
A<I> = e(dQ/de)(2ir/rt) ~ r t ^c{d^l / dr) (2tt /Vt) in angle away from the satellite 
in one orbit (where Q is the angular velocity of the satellite), so a star with 
separation A^ from the satellite, was most likely lost A" or b orbits ago, where 

A* n dR Atfl 

AU ~ ~Z ~ -z • (4) 

27T Ttide d\l ZlT s 

To test how the length of the tidal streamer sample might affect the accuracy of 
the recovered potential we observed 10 random samples of 128 particles in each 
simulation, with the samples chosen to be within A^ = 10, 20, 45, 90, 135 and 
180 degrees of the satellite. The right hand panel of Figure 1 plots the results 
for each simulation as a function of iV or b defined from A\£ via equation (4). The 
plot suggests that if we find ~ 100 stars in a streamer we would like to observe 
them at sufficiently large angle from the satellite to be sensitive to 1-2 orbits. 

Knowledge of Satellite Distance We anticipate knowing the line-of-sight 
velocity and proper motion of the satellite with a precision at least comparable 
to those of the debris stars. However, the satellite's distance r sat may be poorly 
known (or unknown). This means that though we can expect to place strong 
limits on the mass distribution in the Galaxy, Mq^t /r sa ,t) , the physical scale of 
that distribution will be dependent on our measurement of r sat and uncertain 
at the same level. 

Dynamical Friction Binney & Tremaine (1987) show that the frictional 
force on a satellite on a circular orbit in an isothermal sphere is given by (see 
their equation [7.23]) F = —0A28lnAGm 2 Sut /R 2 , where ZnA is the Coulomb Log- 
arithm. The change in orbital energy per orbit will be AEf Yh = F(27ri?)/m sat , 
so 

A F^ Cm 1 

9^ = -0.428/nA^-^27ri?^ — w -2.69/nAs 2 (5) 

e R 2 v 2 irc s 
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Name R dist m sat s Av A/j, A* — 

(kpc) mod M© (km/s) u,as/yr (degrees) 

Sgr 16.0 16.9 5.0E+08 0.151 30.3 252 54 0.393" 

Sex 86.0 19.7 2.6E+07 0.032 6.5 15 11 0.026 

Car 86.6 19.7 1.1E+07 0.024 4.8 11 9 0.016 

LeoII 210.0 21.6 1.4E+07 0.020 3.9 4 8 0.010 



Table 1. Properties of tidal tail data sets associated with each of 
Milky Way's satellites needed to apply the JZSH algorithm. 



Evolution of the Galactic Potential The properties of debris trails will 
also be affected by the time-dependence of the Milky Way's potential (e.g., Zhao 
et al. 1999). However, the existence of the > 7Gyear old thin disk suggests that 
potential evolution has not been significant within the the last iV or b = 1 orbit of 
any of the Galactic satellites. 

4. Conclusion 

Table 1 lists the required accuracy of line-of-sight velocity and proper motion 
measurements of tidal stream stars (columns 6 and 7), length of tidal stream 
population (column 8) and effect of dynamical friction (column 9) for a sample 
of Galactic satellites. The numbers demonstrate that the necessary measure- 
ment accuracies are within the capabilities of SIM and GAIA and that samples 
extending tens of degrees along a streamer to select targets will be sufficient. 
Dynamical friction needs only to be modeled carefully for the innermost satel- 
lites. We conclude that the JZSH algorithm, combined with data from SIM or 
GAIA will significantly improve our knowledge of the Galactic potential. 
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